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SOLUTION  TECHNIQUES  FOR  THE  TRAFFIC  ASSIGNMENT  PROBLEM 


by 

Stephen  P.  Bradley 


1.  INTRODUCTION 

The  object -of  -this  paper  is— to- Investigate  the  application  of  some  large 
scale  solution  techniques  to  the  so  called  traffic  assignment  problem,  or 
multicommodity  flow  problem.  As  this  is  a  network  problem,  there  are  two 
equivalent  formulations.  The  first  is  based  on  the  node-arc  incidence  matrix 
and  is  the  usual  approach.  The  second  is  based  on  an  arc-chain  matrix  and  the 
formulation  is  an  extension  of  Ford  &  Fulkerson's  multi-commodity  work 
Decomposition  techniques  wiUr-be  applied  to  the  first  formulation  and  a  minimum 
chain  network  approach  wtH~be  applied  to  the  second.  Then  an  extension  of  the 
arc-chain  formulation  wtirl  show  it  is  equivalent  to  decomposition  theory  being 
applied  to  a  slightly  different  node-arc  formulation.  , 

The  great  difficulty  in  handling  the  traffic  assignment  problem  is  the 
tremendous  number  of  equations.  If  we  are  to  consider  a  network  with  1,000 
arcs,  300  nodes,  and  100  centroids  (or  commodities),  where  the  centroids  are  a 
subset  of  the  nodes  through  which  all  flow  enters  or  leaves  the  network,  we  are 
faced  with  a  linear  program  with  31,000  equations  and  101,000  variables.  Thus, 
techniques  which  are  not  prlnciply  network  oriented  must  be  discarded.  Both 
compact  triangularizatlon  [2]  and  generalized  upper  bounding  [3]  have  been 
investigated  and  found  wanting  for  the  general  problem,  although  these 
techniques  may  still  prove  to  be  useful  in  solving  the  restricted  master. 
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2.  NODE-ARC  FORMULATION 


The  node-arc  formulation  of  the  traffic  assignment  problem  has  been 
extensively  explained  by  Niels  Jorgensen  [6].  Therefore,  It  will  be  assumed 
that  the  motivation  for  the  formulation  is  understood  so  that  it  will  only  be 
necessary  to  treat  the  problem  from  the  mathematical  standpoint.  The  statement 
of  the  problem  Is  as  follows: 
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The  k  commodity  Is  defined  as  the  flow  that  originates  at  the  k  node.  Thus, 

(10  00  th 

Efv  '  -  q  1  represents  in  matrix  form  the  conservation  equations  for  the  k 

(k) 

commodity  In  terms  of  E  ,  the  node-arc  incidence  matrix,  f^  ,  the  flow  of  the 
kth  commodity  on  arc  ,  and  q ^  ,  the  exogenous  flow  of  the  k**1  commodity  at 
node  N^  .  Note  that  q£^  >  0  q^  <  o  ,  for  1  i  k  .  That  is  the  k^  commodity 
flows  Into  the  network  at  node  Nfc  and  out  of  the  network  at  a  specified  subset  of 
the  remaining  nodes.  The  last  group  of  equations  of  (2.1)  are  the  capacity 

w  00 

constraints,  ie.  \  f.  +  s.  «  k.  .  This  means  that  the  total  flow  on  arc 

v  1  11 

r  (k)  * 

A.  ,  f .  ■  l  f '  ,  must  be  less  than  or  equal  to  the  capacity  of  the  arc  A.  . 

1  1  k  1  1 

The  objective  is  to  minimize  total  travel  time  on  the  network,  thus 

I  t.f.  or  min  Z  —  Y  7  t  f^  . 

i  J J  ik  1J 


min  Z 


Note  that 


>  0  ,  8i  >  0 


Vi,  k 


The  technique  to  be  employed  is  the  general  decomposition  theory  of  Dantzig  [1] 


with  network  techniques  utilized  on  the  subproblems.  For  simplicity  only  two 
subproblems  will  be  considered,  but  the  method  obviously  generalizes  to  any 
number  of  subproblems.  Our  problem  with  only  two  subproblems  is  thus: 
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find  f(1) 

>  o  , 

f*2*  >  0 

,  s  >  0  , 

subject  to 

tf(D 
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If  we  now  consider  the  system  of  equations  Ef  ■  q  ,  it  is  well  known  that 
any  feasible  solution  of  che  system  may  be  represented  by  an  appropriate  weighting 
of  all  the  extreme  point  solutions  of  the  system.  Thus,  we  have  for  the  two 
subproblem** 
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Hence,  letting  f  -  -Z  ,  any  solution  of  (2.1*)  may  be  represented  as  follows: 
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Note  that  any  solution  A^  ,  A^  to  (2.3)  determines  an  f^  ,  f^  by 
equations  (2.2).  We  denote  the  indicated  linear  transformations  as 

Si1)  “  AfiX)  ;  Si2>  “  Afi2)  (2.4) 
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Thus,  we  have  a  problem  that  la  equivalent  to  (2.1*). 


find  A^  >  0  ,  A^^  >  0  and  max  f 
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The  multipliers  associated  with  each  equation  are  Indicated  at  the  left. 

Note  that  beginning  with  equation  (2.3),  the  usual  notation  for  the  decomposition 
approach  has  been  Introduced,  fo  is  simply  equal  to  -Z  so  that  we  maximize 
fo  Instead  of  minimizing  Z  .  The  indicated  matrices  A  and  B  are  simply: 
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The  program  (2.5)  Is  refered  to  as  the  full  master  because  it  contains  all  the 
extreme  point  solutions  of  the  subproblems.  We  will  refer  to  a  restricted  master 
which  is  a  subset  of  the  full  master  that  includes  the  current  basis. 

Assume  for  the  moment  that  we  have  a  basic  feasible  solution  to  (2.5).  That 
Is,  we  have  a  subset  of  k  columns,  of  the  K  columns  of  subproblem  one,  and  a 
subset  of  l  columns,  of  the  L  columns  of  subproblem  two.  Together  with  the 
fo  column,  these  k  +  t  columns  constitute  a  basis  for  the  full  master.  If  we 
then  associate  multipliers  with  each  equation,  we  may  determine  their  values  in 


I 
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the  usual  manner: 
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0  0  (1)  (1)  0  ( 
*  fo  -  1  ;  TT  s}1'  -  y$ L)  ;  tt  sj 


(2) 
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wo 


o  (1)  (2) 

Now  that  we  have  the  multipliers  tt  ,  -wq  ,  -wo  we  must  check  to  see  if  we 
have  an  optimal  solution.  That  is,  we  want  to  examine  the  columns  of  (2.5)  that 
have  the  lowest  cost.  We  must  determine 
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But  we  see  that 


0  (k) 
min  tt  S. 
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0  (k) 
min  y  fK ' 


subject  to 

Ef(k)  -  q*k*  ;  f(k)  >  0 


which  implies  we  must  solve  the  subproblems  of  the  form 

find  f(k)  >  0  ,  min  Z(k) 

subject  to  y  f^  -  (2.6) 

Ef(k>  -  q(k> 

0  0 

where  y  ■  tt  A 

(k)  (k) 

If  the  min  L  <  wo  »  its  relative  cost  factor  is  less  than  zero  and  is 

therefore  a  candidate  to  enter  the  basis.  Since  there  are  in  general  M  sub- 

problems,  we  should  introduce  a  new  column  SA  according  to  the  normal  simplex 

0  (k)  (k}  (k)  (k) 

rules,  ie.  S*  corresponding  to  min  [tt  S*  -  wo  '  ]  •  If  min  *  ■  wo 

k 

for  all  k  ,  we  have  the  optimal  solution. 

Consider  now  the  subproblem  that  must  be  solved. 
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Thus,  our  subproblems  at  each  Iteration  are  the  same  with  a  modified  objective 
function.  The  subproblem3  are  Indeed  merely  minimal  cost  flow  problems  without 
any  capacity  constraints.  These  may  easily  be  solved  by  network  techniques.  A 
straight-forward  application  of  the  out-of-kilter  algorithm  will  yield  a  ready 
solution  and  has  the  advantage  that  the  solution  to  the  previous  subproblem  may 
be  used  as  a  starting  solution  to  the  next,  thus  speeding  convergence.  Another 
solution  technique  is  the  dynamic  programming  algorithm  for  the  minimal  spanning 
tree  routed  at  the  k^  node.  With  this  method,  it  is  also  necessary  to  allocate 
flow  according  to  the  demands  at  each  node  and  sum  the  resulting  flows  ..or  each  arc. 
Thus,  we  see  that  we  have  in  general  M  uncapacitated  minimal  cost  flow  problems 
to  solve  in  conjunction  with  a  restricted  master. 

A  few  comments  on  the  interpretation  of  the  objective  function  are  in  order. 

0 

Since  itq  -  1  ,  the  objective  of  the  subproblems  is 
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Y  f 
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Therefore,  each  arc  has  a  "cost"  associated  with  it  which  is  the  sum  of  the  given 

0 

unconstrained  travel  time  t?  ,  and  an  additional  tarrif  tt^  due  to  the  demand 
for  the  arc  exceeding  its  capacity.  In  effect,  we  put  a  price  on  the  arcs  with 
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insufficient  capacity  to  discourage  flow  from  using  these  arcs.  If  the  arc  is  not 
at  capacity,  we  have  more  capacity  than  is  demanded,  and  thus  the  price  on  the  arc 
is  zero.  At  each  iteration  the  prices  are  re-evaluated  and  flow  rerouted  accordingly 
until  the  optimal  is  reached.  Note  that  the  objective  functions  are  the  same  for 
all  the  subproblems.  This  is  consistent  with  not  discriminating  between 
commodities  on  any  arc.  Thus,  the  price  fo1*  any  commodity  to  use  a  particular  arc 
is  the  same  for  each. 

The  procedure  is  then  to  introduce  the  winning  column  into  the  restricted 
master,  which  consists  of  the  current  basis  and  the  nonbaslc  slacks.  We  then 
optimize  the  restricted  master,  yielding  a  new  basis.  A  new  set  of  multipliers 
may  then  be  calculated  and  the  subproblems  resolved  with  these  new  prices,  etc. 

Since  at  each  stage  we  are  in  effect  optimizing  the  restricted  master,  the  prices 

0 

associated  with  the  nonbaslc  sleek  variables  will  be  >  0  .  Note  the  prices 

0 

associated  with  basic  slack  variables  are  -  0  . 

The  solution  of  the  subproblems,  although  time  consuming  due  to  the  large 

number,  does  not  seem  to  be  the  restrictive  aspect  in  this  approach.  Since  these 

problems  must  be  solved  even  in  the  uncapacitated  traffic  assignment  problem,  it  is 

not  surprising  to  find  they  must  be  solved  in  the  capacitated  case.  At  each 

iteration,  ie.  solution  of  the  subproblems,  we  must  find  the  new  prices  to  add  on  to 

the  arcs  at  capacity.  Since  we  know  that  when  alack  occurs  on  arc  A^  the 

0 

corresponding  price  ir^  is  equal  to  zero,  we  need  only  be  concerned  with  the 

(k) 

prices  associated  with  the  arcs  at  capacity  and  the  multipliers  -po  associated 
with  the  subsum  equations.  Thus,  at  least  theoretically  it  is  only  necessary  to 
Invert  a  matrix  of  this  rank.  As  the  inversion  process  limits  the  size  of  the 
network  problem  we  can  handle,  we  would  like  to  find  avi  alternative  method  for 
finding  the  prices. 


-as f 
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3.  ARC-CHAIN  FORMULATION 

Instead  of  dealing  with  the  usual  node-arc  formulation,  it  is  well  known 
that  there  is  an  equivalent  formulation  in  terms  of  an  arc-chain  matrix.  Since 
an  interesting  method  of  solution  is  available  in  this  formulation,  it  will  be 
presented  here.  First  it  is  necessary  to  introduce  some  concepts  relevant  to  the 
arc-chain  formulation.  Let  A  ■  |i..|  be  an  arc-chain  matrix.  Then  we  define 


M 


‘ij 


1  if  arc  is  in  chain  Cj 
0  otherwise 


Now  for  the  traffic  assignment  problem,  we  will  change  the  definition  of  a 

commodity  slightly.  A  commodity  is  defined  to  be  flow  between  a  particular  origln- 

(k) 


raooi 

_aij  j 


be  the  arc-chain  matrix  consisting  of  all 


destination  pair.  Let  A 

chains  connecting  the  kth  origin-destination  pair.  In  the  previous  formulation 

(k) 

we  had  flow  on  an  arc,  here  we  have  flow  on  a  chain.  Let  be  the  flow  of 

the  kth  commodity  on  chain  Cj  .  We  may  now  represent  the  capacity  constraint  on 
the  1th  arc  as  follows: 


(l)v(D  .  V  _ (2)v(2) 
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+  ...  +  l  a^'xr'  +  s4  -  K, 


(3.1) 
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i  m  1,2,  . . . ,  n 

o 

This  merely  means  that  the  flow  in  arc  A^  is  equal  to  the  sum  of  the  chain  flows 
through  the  arc  A^  ,  and  this  total  flow  plus  the  slack  must  equal  the  capacity 
of  the  arc  A^  . 

The  second  requirement  of  the  problem  is  that  there  is  a  specified  amount  of 
flow  between  any  origin-destination  pair.  Therefore,  the  total  chain  flows 
associated  with  the  kC^  origin-destination  pair  must  equal  the  required  flow. 


l  x.(k)  -  v(k) 

Jk  Jk 


Vk 


(3.2) 
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(k) 

Of  course,  all  Che  chain  flows  must;  be  nonnegative,  Xj  >0  . 

(k) 

Finally,  we  define  c^  as  the  total  travel  time  on  chain  Cj  between 
the  kfck  origin-destination  pair.  Since  we  wish  to  minimire  total  travel  time 
on  the  network,  our  objective  function  is 


min  Z 


l  cj^xj15  +  l  cj2)Xj2)  +  ....  +  l  cJ(o)Xj(n) 

h  h  Jm 


(3.3) 


Thus,  we  have  formulated  the  traffic  assignment  as  the  following  linear  program. 
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An  extension  of  the  multi-commodity  work  of  Ford  &  Fulkerson  [5]  leads  to 
a  solution  of  the  problem.  Assume  for  the  moment  we  have  a  known  basic  feasible 

0  (k ) 

solution  to  (3.4).  Thus,  we  can  calculate  ,  -pq  as  we  did  in  the  previous 
section.  Note  that  these  prices  are  the  same  type  as  those  used  in  the  decomposition 
approach,  thus  their  sign  is  opposite  that  of  the  usual  simplex  multiplier.  Let  us 
now  consider  whether  the  prices  are  optimal  for  the  program.  Pricing  out  we  have 


*1 


A(k) 

i  ij 


>  0 


VJ  ,  k 


(3.5) 
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for  the  optimality  criterion.  Notice  that  to  test  optimality  It  is  not 


-(k) 

necessary  to  look  at  all  c  ^ 


_(k) 

If  we  look  at  c  *  min  c'  ,  we  know  if  we 

8  M_  J 


are  optimal  by  whether  c  is  nonnegative  or  not;  if  c  >  0  ,  we  are  optimal 

ft  8  ® 


Also  note  that  cj^  t  the  total  travel  time  on  chain  between  the  kth  origin- 

destination  pair  may  be  represented  as  follows: 


00 


t(k) 

ij 


where  t^  is  the  travel  time  on  arc  just  as  before.  The  total  travel  time 

on  a  chain  equals  the  sum  of  the  travel  times  on  the  arcs  in  that  chain.  Thus, 
we  may  write 


(k) 
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ij 


0  0 

t^  >  0  and  we  will  show  >  0  ,  thus  (t^  +  tt^)  >  0  .  If  we  now  consider 
0 

(Ci  +  88  t*ie  °i  fcbe  arcs,  we  may  interpret  the  minimum  relative  cost 

factor  in  an  interesting  way. 
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Since  uo  i8  a  constant  for  a  particular  origin-destination  pair,  we  need  only 
look  at 


min  l  (t  +  ir.jaj^ 
j  i  1 
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We  see  that  this  merely  represents  the  shortest  distance  between  the  k^  origin- 

0 

destination  pair,  when  (t^  +  is  considered  as  the  length  of  the  arc. 

Therefore  to  check  optimality,  we  need  only  find  the  shortest  paths  associated  with 

(k) 

the  origin-destination  pairs,  subtract  the  appropriate  po  and  take  the  minimum 
over  all  k  . 


c 
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min  •  min  £(t  +  n.)a^l  -pg^ 

k  j  [i  2  J 


(3.8) 


Thus,  we  see  that  is  is  not  necessary  to  explicitly  enumerate  all  the  chains  of  the 
network.  Only  the  shortest  chain  for  each  commodity  must  be  examined.  Not  having 
to  know  the  arc-chain  matrix  explicitly  is  absolutely  necessary  for  this  method  to 
be  computationally  feasible. 

A  few  comments  on  the  signs  of  the  multipliers  are  now  in  order.  We  have 
defined  the  multipliers  consistent  with  the  decomposition  theory,  and  thus  they 
are  opposite  in  sign  with  respect  to  the  usual  simplex  multipliers.  The  multipliers 
associated  with  the  current  basis  of  the  restricted  master  are  defined  by 

0 

tt  Pj«0  j/0;  tt'Pq  ■  1  for  basic  columns 

0  0 

It  is  easy  to  see  that  if  s^  is  basic  *  0  .  Now  if  s^  is  nonbaslc  >  0 
since  we  are  at  each  stage  optimal  with  respect  to  a  subset  of  columns  including  the 

current  basis  and  the  nonbaslc  slack.  Note  this  is  exactly  analogous  to  the 

0 

decomposition  situation.  The  are  in  fact  prices,  just  as  they  were  before, 

and  represent  additional  tarrifs  associated  with  the  use  of  the  capaclted  arcs. 

As  before,  this  additional  tarrif  is  added  to  the  unconstrained  travel  time  t^ 

forming  a  new  "cost",  and  the  flow  is  rerouted  accordingly.  We  may  also  determine 
(k) 

the  sign  of  -uq  in  an  easy  manner.  Since  at  least  one  variable  from  each 
commodity  subset  must  be  basic  [3],  we  may  write 

-PO^  +  £(t^  +  n^Ja^  ■  0  for  some  j 


( 

j 

i: 

t 

i 


) 


**  wa  j  *  M 
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0 

or  rewriting  and  noting  that  t^  >  0  and  tt^  >  0  ,  we  have 

-uok)  "  I(tt  +  >0  Vk  (3.9) 

UO^  is  the  length  of  the  shortest  path  between  the  kC^  origin-destination  pair 
that  is  currently  in  the  basis.  All  such  paths  for  the  kC^  commodity  must  have  the 
same  length. 

Rather  than  discussing  the  relative  merits  of  this  arc-chain  formulation  at 
this  point,  we  will  put  this  aside  until  the  end  of  the  next  section.  In  the 
next  section,  we  will  show  that  the  arc-chain  formulation  is  equivalent  to  a 
decomposition  formulation  and  thus  the  problems  encountered  will  be  similar. 
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4.  EXTENSION  OF  THE  ARC-CHAIN  FORMULATION 

The  similarity  between  the  two  previous  formulations  lead  us  to 
investigate  further  their  relationship.  The  outstanding  similarity  is  that  in 
each  case  prices  are  put  on  the  arcs  at  capacity  which  are  added  directly  to  the 
uncapacitated  travel  times  of  the  arcs.  The  flow  allocation  is  then  re-examined 
with  respect  to  these  new  composite  travel  times.  Also  in  the  decomposition 
approach  we  solve  a  subproblem  which  is  a  minimal  spanning  tree,  or  shortest 
route  to  all  nodes  of  the  network.  In  the  arc-chain  formulation,  we  solve  a 
shortest  route  between  an  origin-destination  pair.  Let  us  then  make  the 
following  transformation  of  variables  on  the  arc-chain  formulation  (3.4): 


Y(k)  _  .(k).,(k) 

XJ  V 


Vk 


(4.1) 


This  results  from  dividing  each  flow  requirement  equation  by  the  required  flow; 
in  effect  normalizing  the  equations.  Thus,  we  may  write  the  arc-chain  problem 
as  follows  in  terms  of  the  new  variables. 


min 


I  C.(1)A.(1)  +  V(2)  l  C<2\<2>  +  ....  +  V(m)  J  cJm;A(m) 
JlJ  J  ^  J  *mJ  * 


(m)  r  _(m), 


subject  to 


V(1)  I  +  V(2)  l  a(2M2)  + - +  V(m)  l  aj”V(m)  + 

J*J  J  4  J  4  *J  J 

1  J2  3ro 


s .  *  k. 


Vi 


l  a!1’ 
h 


-  1 


(4.2) 


l  a.<2) 

J, 3 


■  1 


l  A,0'0 
V 


-  1 


MW1  W'V *>  _ 
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Now  our  program  looks  exactly  like  the  restricted  master  of  a  decomposition 

(k)  (k) 

program.  Note  that  each  subsum  of  A^  equals  one  and  thus  the  Aj  mny 
be  considered  as  weights  associated  with  the  corresponding  columns.  Also,  it 
is  easy  to  see  that  each  column  is  a  path  between  the  k^  origin-destination 
pair  with  the  appropriate  flow  on  it,  namely  V  .  Thus,  equation  set  (4.2) 
is  really  a  full  master  for  a  decomposition  problem  whose  subproblems  are  shortest 
route  problems  between  an  origin-destination  pair.  By  full  master,  we  mean  that 
(4.2)  contains  all  the  extreme  point  solutions  of  the  subproblems.  The  reason 
why  we  do  not  have  to  explicitly  know  the  arc-chain  matrix  is  simply  that  we 
bring  up  the  needed  extreme  point  solution,  ie.  shortest  path,  from  the  sub¬ 
problem  only  when  it  prices  out  negative.  Since  network  techniques  are  used  on 
the  subproblems,  they  could  be  formulated  in  the  usual  node-arc  manner.  Thus, 
we  have  shown  that  the  arc-chain  formulation  and  solution  technique  is  equivalent 
to  the  decomposition  technique.  The  decomposition  problem  for  which  (4.2)  is  the 
full  master  is  as  follows: 

(k) 

find  i  >0  s  >  0  ,  and  min  Z 


subject  to 


t  t  t  t 
E 


I  I  I 


1 03 


E 

I  I 


Z  (min) 


(1)" 

rq(i)i 

(2) 

q(2> 

(3) 

q<3> 

(P) 

s 

,<P> 

K 

-  J 

where  q 


(k) 


-v 

0 


.00 


(k) 


(4.3) 


The  similarity  between  this  and  our  original  node-arc  formulation  is  very  great. 
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Where  before  we  had  M  usbproblems ,  which  involved  finding  the  shortest  path  to 
all  other  nodes,  here  we  have  M(M-l)  subproblems,  which  involve  finding  the 
shortest  path  between  an  origin-destination  pair. 

Thus,  we  have  shown  that  the  arc-chain  formulation  simplifies  the  subproblem 
considerably,  but  at  the  same  time  it  expands  the  number  of  subproblems 
considerably.  Both  the  node-arc  and  the  arc-chain  procedures  depend  on  the  same 
principles.  We  solve  a  number  of  routing  problems  given  a  set  of  prices  on  the 
capacited  arcs  and  check  optimality.  If  we  are  not  optimal,  we  re-evaluate  our 
prices  and  resolve  our  routing  problems  again  with  respect  to  the  new  prices.  In 
both  procedures  the  prices  associated  with  the  arcs  not  at  capacity  are  zero. 

Thus,  in  the  aic-chain  formulation  we  must  theoretically  invert  a  matrix  the 
rank  of  the  number  of  arcs  at  capacity  plus  the  number  of  requirement  equations. 

A  criticism  that  might  be  lodged  against  the  arc-chain  formulation  is  that  the 
number  of  subproblems  has  become  prohibitively  large.  Indeed  the  number  of 
subproblems  places  a  bound  on  the  size  problem  we  can  handle,  but  in  fact  these 
M(M-l)  shortest  route  problems  may  be  solved  more  quickly  than  the  M  un¬ 
constrained  minimal  cost  flow  problems.  This  is  true  since  the  M(M-l)  problems 
are  solved  by  applying  the  dynamic  programming  algorithm  for  minimal  spanning  tree 
routed  at  a  particular  node  M  times.  Each  application  solves  M-l  of  the  sub- 
problems.  Since  this  must  be  done  in  the  minimal  cost  flow  situation  plus 
assigning  flows  and  adding  appropriately  for  each  arc,  it  will  be  quicker.  The 
difficulty  with  the  method  is  that  there  are  considerably  more  multipliers 
associated  with  subsum  equations  to  be  calculated. 

Note  that  in  each  case  we  assumed  that  we  had  a  basic  feasible  solution  to 
the  restricted  master.  This  is  usually  not  the  case,  but  the  method  outlined  by 
Dantzlg  [1]  may  be  used  to  obtain  one.  We  simply  add  artificial  variables  to  the 
restricted  master  and  minimize  the  infeasibility  form.  The  method  is  the  same  as 
that  described  for  the  decomposition  approach.  Here  the  unconstrained  travel  time 
on  the  arcs  is  considered  zero,  and  a  price  is  put  on  those  arcs  that  are  associated 


-—Trv&ssssi 


mrvsssr . 


with  an  infeasible  equation. 

fius,  we  have  shown  how  a  pricing  mechanism  forces  the  re-allocation  of 
commodity  flow  until  an  overall  optimum  is  attained.  For  a  given  set  of  prices, 
we  check  to  see  if  the  current  flow  allocation  is  optimal;  and  if  not  we  find 
new  prices  for  the  arcs  and  reroute  flow  accordingly.  The  simplex  criteria  tell 
us  which  new  route  (chain)  should  be  considered  and  which  old  route  should  be 
dropped.  Then  new  prices  are  found  consistent  with  the  new  basis  and  the 
processes  repeats  itself. 
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5.  COMPUTATIONAL  STATUS 

Now  that  the  theory  has  been  demonstrated,  it  seems  appropriate  to  say  a 
few  more  words  concerning  computational  possibilities.  We  have  indeed  shown  that 
the  problem  can  be  formulated  as  two  slightly  different  decomposition  problems. 

The  first  defines  commodity  as  that  flow  leaving  the  kC^  node.  Thus,  for  M 
commodities,  or  centroids,  we  have  M  uncontrained  minimal  cost  flow  problems. 

The  second,  an  extension  of  Ford  &  Fulkerson's  work,  defines  commodity  as  flow 
between  the  kC^  origin-destination  pair.  Thus,  for  M  centroids  we  have  M(M-l) 
shortest  route  problems.  The  differences  are  only  slight  and  only  slightly 
different  solution  techniques  are  implied. 

The  unconstrained  minimal  cost  flow  problems  of  the  first  formulation  may 

be  solved  in  two  ways.  First,  a  straight  forward  application  of  the  out-of-kilter 

algorithm  will  yield  a  rapid  solution  as  the  problems  are  unconstrained.  Also, 

since  the  objective  function  is  the  same  for  each  subproblem  at  any  iteration, 

th 

the  solution  of  the  i  subproblem  can  be  used  as  a  starting  solution  for  the 
i  +  1^  subproblem  at  any  iteration  -  thus  speeding  convergence.  As  an 
alternative,  it  appears  that  using  the  solution  of  the  I1*1  subproblem  at 
iteration  t  as  the  starting  solution  for  the  i^  subproblem  at  iteration 
t  +  1  will  make  convergence  extremely  rapid  as  the  right  hand  sides  will  be 
identical  and  the  objective  functions  only  slightly  modified.  Of  course,  this 
entails  storing  in  some  manner  the  M  solution  vectors  of  the  subproblems  from 
Iteration  to  Iteration.  Note  that  this  causes  no  real  difficulty  as  the  solution 
vectors  currently  in  the  basis  of  the  restricted  master  will  be,  if  not 
equivalent,  just  as  useful.  The  second  method  to  solve  the  subproblems,  which 
has  been  used  in  practice,  is  to  solve  a  minimum  tree  routed  at  a  particular  node 
for  each  subproblem.  This  is  done  extremely  rapidly  by  application  of  the  dynamic 
programming  algorithm  for  "shortest  path"  (actually  shortest  path  to  all  nodes 
and  thus  minimum  tree).  This  computation  may  be  shortened  by  judicious  choice 
of  the  order  in  which  the  subproblems  are  solved.  That  i9,  if  the  appropriate 


node  is  chosen  adjacent  to  the  one  for  which  the  minimum  tree  was  just  calculated, 

a  great  deal  of  the  new  minimum  tree  will  be  the  same.  Unfortunately,  their 

presently  exists  no  systematic  technique  for  doing  this.  To  find  the  column 

to  introduce  into  the  basis,  we  need  the  flows  on  the  arcs.  Thus,  for  each 

destination  node  of  a  subproblem,  we  traverse  the  shortest  route  to  the  origin 

00 

(ie.  in  the  reverse  direction)  assigning  the  demand  on  each  arc  . 

When  this  is  done  for  each  destination  node,  we  simply  sum  the  resulting  demands 
on  each  arc  to  obtain  the  required  total  flow  on  that  arc.  Although  this  method 
of  solution  is  more  difficult  to  describe  than  the  out-of-kilter  method,  many 
argue  that  it  is  as  efficient. 

Now  with  the  subproblems  solved  at  some  iteration,  we  want  to  introduce 
the  appropriate  column  into  the  restricted  master.  We  could  use  the  straight 
forward  application  of  a  linear  programming  code  to  optimize  the  restricted 
master,  as  this  yields  the  desired  prices  directly.  Clearly,  if  this  method  is 
used  all  those  columns  that  currently  price  out  negatively  from  solution  of  the 
subprobleras  should  be  introduced  into  the  restricted  master;  then  the  restricted 
master  is  optimized.  The  number  of  equations  is  then  equal  to  the  number  of  arcs 
plus  the  number  of  commodities  (subproblems).  Thus,  the  rank  of  the  restricted 
master  constrains  the  size  problem  you  can  solve  (ie.  as  large  as  the  linear 
programming  code  employed).  There  is  an  alternative  approach,  although  it  has 
not  yet  been  fully  worked  out.  That  is,  we  utilize  the  fact  we  know  the  price 
is  zero  on  all  arcs  not  at  capacity.  Thus,  if  we  select  our  matrix  to  invert, 
it  must  only  be  of  rank  equal  to  the  number  of  capacitated  arcs  plus  the  number 
of  commodities.  This  is  a  sizable  reduction  as  the  number  of  capacitated  arcs  is 
rarely  more  than  sixty  percent  of  the  total  and  often  as  low  as  thirty  percent. 

This  approach  has  the  obvious  difficulty  that  the  number  of  capacitated  arcs 
may  change  from  iteration  to  iteration. 
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The  second  formulation  is  of  interest  because  the  technique  for  the  sub¬ 
problems  is  straight  forward.  Here  we  have  M(M-l)  shortest  route  problems 
between  origin-destination  pairs.  We  simply  apply  the  dynamic  programming 
algorithm  M  times;  each  application  yielding  the  solutions  to  M-l  of  the 
shortest  route  subproblems.  Note  that  no  back  search  is  necessary  to  find  the 
flows.  The  only  difficulty  with  this  solution  approach  in  general  is  that  the 
restricted  master  has  n  equations  for  the  capacity  constraints  as  before  but 
also  M(M-l)  commodity  equations  instead  of  M.  Therefore,  the  number  of 
commodity  equations  may  become  prohibitive  unless  special  techniques  are 
employed.  For  a  large  problem,  it  is  suggested  that  the  generalized  upper 
bounding  algorithm  of  Dantzlg  and  Van  Slyke  be  utilized  [3].  The  algorithm 
given  in  Vol.  II  of  their  paper  has  a  "working  basis"  of  rank  equal  to  the 
number  of  arcs. 

In  conclusion,  it  would  appear  that  a  reasonable  size  problem  can  be  solved 
straight  forwardly  by  either  method.  The  bound  on  this  "reasonable  size"  is 
placed  either  by  the  linear  programming  code  size  utilized  to  solve  the 
restricted  master  or  by  the  computing  time  utilized  in  solution  of  the  subproblems. 
For  present  computers,  a  chain  job  configuration  would  have  to  be  used  on  a  large 
problem;  this  could  be  quite  time  consuming,  but  future  machines  will  eliminate  this 
difficulty. 

At  this  stage  there  are  three  areas  where  further  research  seems  justified. 
First,  an  approach  using  only  network  techniques  would  be  desirable,  but  the  simple 
pricing  mechanism  Implied  by  the  decomposition  approach  would  seem  hard  to  improve 
upon.  Second,  an  efficient  technique  to  calculate  the  prices  that  uses  only  the 
Information  associated  with  the  capacitated  arcs  would  greatly  Increase  the  size 
problem  the  above  techniques  could  accomodate.  Third,  in  practice  a  systematic 
technique  of  guessing  the  prices  might  be  the  best.  Finally,  as  computer  technology 
Increases  the  amount  of  rapid  access  storage  available  to  the  user,  the  above 
techniques  should  prove  to  be  extremely  rapid. 


APPENDIX 


An  intersting  extension  of  the  classical  traffic  assignment  problem  has 
been  proposed  by  A.  Ridley. +  The  basic  concern  is  to  add  a  budget  which  can  be 

used  to  Increase  the  capacity  of  any  arc.  Thus,  the  capacity  constraints  for  the 
problem  take  on  a  different  form.  The  right  hand  side  now  becomes  the  sum  of  two 
capacities,  the  present  existing  capacity  plus  that  capacity  which  can  be  added 
by  spending  money  on  that  particular  arc.  Thus,  we  have 


(k) 


KJ  +  kj 


on 


k  2 


(k) 


+  s 


Kj+kj 


VJ 


(A.  1) 


Since  we  want  to  minimize  total  travel  time  on  the  network,  we  could  Increase  the 
capacities  on  all  the  shortest  paths;  but  this  is  not  possible  with  only  a  finite 
amount  of  money.  If  we  call  B  ,  the  budget  or  the  total  amount  we  can  spend;  and 
call  r^  the  cost  per  unit  increase  in  capacity  on  the  ith  arc  we  have 

£  r.k  <  B  on  £  r.k  +  s'  (A. 2) 

j  J  J  "  j  J  3 


That  is  the  amount  of  money  spent  on  each  arc,  summed  over  all  arcs,  must  be  less 
than  or  equal  the  total  budget.  The  new  problem  in  matrix  notation  is  then 
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The  reason  for  considering  this  variation  is  to  show  that  the  proposed  solution 
technique  also  solves  this  problem  as  an  easy  extension.  The  subproblems  are 
obviously  exactly  the  same  as  they  were  in  the  initial  formulation.  Thus,  let  us 
now  consider  the  restricted  master. 


find  A^  >  0  A^^  >  0  s>0  k>0  s'  >  0  max  f. 


Vo  +  I  Si1>Ai1)  +  l  S{2)^2)  +  Is  -  Ik  -  K 


(A. 4) 


r  k  +  8 '  ■  B 

-  1 


-  1 


Note  that  this  only  adds  one  more  equation  to  the  problem  whll'  adding  several 
more  variables.  Thus,  any  straight  forward  method  for  optimising  the  restricted 
master,  such  as  the  revised  simplex  method,  will  yield  prices  just  as  it  did 
before.  In  fact,  the  solution  technique  for  this  problem  is  identical  with  that 
of  the  first  formulation  except  shortcuts  for  obtaining  prices  may  not  be 
applicable. 

Since  the  objective  is  to  minimize  total  travel  time  regardless  of  cost, 
assuming  we  are  within  the  budget,  the  solution  is  Intuitive.  We  will  spend 
money  to  increase  the  capacity  of  the  shortest  paths  on  critical  arcs  until  one  of 
two  things  happens.  Either  all  commodities  will  be  traveling  over  their  shortest 
paths  only  or  the  budget  will  be  constraining. 


*.  Vi V 
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